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Abstract 

General formalism of quantum field theory and addition the- 
orem for Bessel functions are applied to derive formula for Casimir- 
Polder energy of interaction between a polarizable particle and a di- 
lute dielectric ball. The equivalence of dipole-dipole interaction and 
Casimir energy for dilute homogeneous dielectrics is shown. A novel 
method is used to derive Casimir energy of a dilute dielectric ball with- 
out divergences in calculations. Physically realistic model of a dilute 
ball is discussed. Different approaches to the calculation of Casimir 
energy of a dielectric ball are reviewed. 
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1 Introduction 



The study of spherical geometry in Casimir effect meets a lot of technical 
problems, the most difficult one is the problem of divergences, which appear 
in many expressions. The special interest to the subject raised after the 
series of articles by Julian Schwinger where he had suggested a connection 
between Casimir effect and sonoluminescence [|T], |2|] . Different regularizations 
were used to obtain finite results in spherical geometry, however, results of 
the calculations strongly depended on the regularization scheme and regu- 
larization parameters (a discussion of divergences and their regularization 
can be found in In the pesent paper we show that it is necessary that 
regularization parameters be present even in final expressions, they have real 
physical meaning and we will clarify it. 

We study a dielectric nonmagnetic ball of radius a and permittivity e, 
surrounded by a vacuum. The ball is dilute, i.e. all final expressions are 
obtained under the assumption e — 1 <C 1. Dispersion is essential in the 
following, so the permittivity e depends on frequency. 

In section 2 we follow the formalism which was developed by E.Lifshitz 
et.al. |4[] and K.Milton et.al. || . We start from a short overview of known facts. 
Then we derive Casimir-Polder type energy of interaction between a dielec- 
tric ball and a particle of constant polarizability a, which is placed at the 
distance r > a from the centre of the ball. The addition theorem for Bessel 
functions is used. Another application of this technique can be found in |J 
where Casimir energy of polarizable particle located in the neighbourhood 
of perfectly conducting wedge was calculated. In section 3 the equivalence 
of dipole-dipole interaction and Casimir energy in the lowest order (e — l) 2 
is proved for homogeneous dielectric media. This fact is used in section 4 to 
calculate non-dispersive contribution to the energy, dispersive contribution 
is briefly discussed. We argue that Casimir surface force is attractive in a 
realistic model of dielectric ball and includes volume and surface terms in 
the order (e — l) 2 . 

Various regularizations are reviewed in the Appendix. 

We put h = c = 1. Heaviside-Lorentz units are used. 
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2 Polarizable particle and a ball 



The change in the ground state energy E of the system under the infinitesimal 
variation of e is 

-+oo 

■ )lr ,1 '< 

-oo 
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5E = 2 J J ^ fe ( x '^) jD pp( £ ' r ' r ' u; )- 



Here Djk(e, r, r, u) is a Fourier component of electric field propagator trace, 
it satisfies equations 



e(r, u) lo 5 jm - rot jZ rot Zm 



D mk (e, r, r', u) = -u 2 5{r - r')6 jk (2) 



The system of equations for this Green's function was discussed extensively 
in 0, |5| (there is a discussion of these equations in the next section). The 
solution of this system for spherical geometry with standard boundary con- 
ditions classically imposed at r = a can be written as in @ (^-functions are 
omitted since we are interested in the limit r — > r'): 



oo I 

D jh (e,v,r',uj) = -Y.H (" 2 Fi(r,r')X jlrn (m* k imW) + 

1=1 m=—l 



+ ^rot r rot r ,G z (r,r / )X jZm (fi)X* im (fi'))- (3) 

Here we have used the following notations (Xu m (Q) are vector spherical har- 
monics; ji(r),fi^\r) are spherical Bessel functions, e z (r) = r/ij 1 ' '(r),s z (r) = 
rji(r) are Riccati-Bessel functions ||): 

X am (Q) = 2_ (LY lm (m (4) 
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differentiation is taken over the whole argument. 

When the point particle of constant polarizability a is inserted into the 
point r, |r| > a from the centre of the ball, the energy change is given 
by (H|) with 5e = 47ra5 3 (r — x). However, the contact terms (empty space 
contribution) have to be subtracted from ([!]), i.e. when physical quantities are 
calculated in the region r,r' > a the term z|uj|/i[ (|aj|r>)jj(|uj|r<) has to be 

subtracted from ([5]) (for r, r' < a in full analogy the term ikji(kr < )h^\kr > ) 
should be subtracted from (||)). Doing so, we have to substitute F h G\ instead 
of Fi, Gi in all expressions, where 



Fi, Gi 



-iApakjiikr^jiikry), k = \u\y/e, r,r' < a, 



iB F ^ G \(j\hl 1 \\uj\r > )h\'' > (luolr^ , r,r' > a 



The Casimir-Polder energy of this configuration is 
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We perform a Euclidean rotation then: u 
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Let x — Loa . For E\ (r, a) we obtain 
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This expression can be transformed to a simple formula in the limit e— 1 <C 1. 
The functions Bp and Bq are proportional to (e— 1) in this limit. To proceed, 
the following addition theorem for Bessel functions M is useful: 



u 



(p, k,x,p) = ^(21 + l)si(xp)ei(xk)Pi(cos( 
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xe xp pk 
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To our knowledge this formula was first used in Casimir effect calculations in 



10 1, where it was applied to analytic calculation of Casimir energy of perfectly 



conducting spherical shell and dilute dielectric ball satisfying efj, — 1. 

In our case it can be applied as follows. The simple identity holds (k > 
p > for definiteness) : 



dx'^^(2l + l)f{x)si{xp)ei{xk)si{xp)ei{xk) 
I f k+p dpp 
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Only the first order ~ {e — 1) is needed in E\. We put k = r /a, p = 1 and 
use ([18]) and its obvious generalizations in (|15D to calculate E\. Finally we 

get 



23 e - 1 a 3 (5r 2 + a 2 ) 
i(r, a)--— a - _ 



r > a . 



(21) 



Substitution e — 1 = 4:TrN mo iab a ii in the limit r ^> a yields 

/47ra 3 \ — 23aa, 
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Thus in this limit the famous Casimir-Polder energy of interaction between 



two polarizable particles E Cas -Poi [0] can be obtained directly from (p2|). 

Certainly, this is not simply a coincidence, and next section is devoted to 
the proof of equivalence of dipole-dipole interaction and Casimir energy in 
the order (e — l) 2 . 



3 Dipole interaction and Casimir energy 

We define differential operator 



L jm (e,r,iuj) 



s r 



u 2 5 jm + rot jiiotim 



(23) 
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Figure 1: dielectric area is shaded; point particle is located at point r' 

Equations for electric field propagator are well known due to works [|]. : 

L jm (e, r, iu)D mk (e, r, r', iu) = -u 2 5(r - r')S jk . (24) 

Suppose that e has changed in the point r 2 . The changed system satisfies 
equation 

L jm (e + 8e(i\u\) 8(r 2 - r),r, iu) D m k{s + 8e(i\u\)8(r 2 - r),r,r',iu) = 

= -u 2 8(r-r')8 jk (25) 

Let us define 8Djk(r, r', iuo) as a term proportional to fe(i|a;|) in the expansion 

D jk = D jk + 5D jk (v, r', iu) + 0((fe(*M)) 2 ) (26) 

If we consider equation ( |25D and write series in powers of constant 5e(i|o;|), 
in the first order we get the equality: 

8e(i\u\)8(r 2 - r) u 2 D jk (e,r, r',iu) + L jm (e, r, iu)8D mk (r, r', iu) = 0. (27) 

We multiply equation (|27|) by D P j(r", r, iu) on the left , integrate over r and 
sum over j. Using properties D p j(r", r, iu) = Dj p (r,r" ,iu) , L^ m = L m j and 
equation (fj4|) it is straightforward to obtain 

8D pk (r", r', iu) = D pj (e, r", r 2 , iu)5e(i\u\)D jk (e, r 2 , r', iu). (28) 
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Now suppose that D m k satisfies (p4|) with e — 1 in all space. Then 
we insert dielectric of arbitrary shape with constant dielectric permittivity 
e(i|u;|) so that 5e(i|u;|) = e: (i | a; | ) — 1. For this case 5-function in ( P7j ) has to 
be substituted by function equal to outside the dielectric and equal to 1 
inside the dielectric. Denoting D^(r", r', iu) = D pk (e = 1, r", r', iu) , we get 
the formula 

D pk (e,r'\r'Ju) - D$(r",r',iu) = 

= J d 3 v 2 D^(v",v 2 ,iu;)^m) - l)Df k \r 2 ,v'M + 0((e(i\u\) - l) 2 ) 

(29) 

Figure 1 illustrates formula fl2"9"|). When r" r' the formula for Casimir 
energy of interaction between a polarizable particle located at point r' and 
dielectric can be written as 



/-t-oo 
du ai(r', i\u\)(D pp (e, r", r, iu) - D p °J (r", r ', iu)) 
-oo 

= J d 3 r 2 p(r 2 )£/(r',r 2 ) + 0(a 3 ), 



(30) 



where 

- + 0O 



/+oo 
duA<Ka x {r', i\u\)a 2 (r 2 , i\u\) (^(r', r 2 , iu)) 2 . (31) 
-oo 



Equation for potential (|3"T| ) can also be obtained from dipole-dipole in- 
teraction of two neutral atoms |12] . We refer the reader to this book. If we 
denote r = |r 2 — r'|, then (see |]l2"fl) 



rCO r 9 ^ fi 



U(r) = - / uj a\(iuj) a 2 (iu) e 



ur (oor) 2 (ur) 3 (cur) 4 

(32) 



du 



If we define atomic absorption frequency as u , then A ~ l/tuo is a char- 
acteristic absorption wavelength in the spectrum of interacting atoms. For 
distances r 3> A the main contribution to ([!!]) gives the frequency region 
uj < 1/r <C 1/A. In this frequency region according to general properties 
of dielectric permittivity on imaginary axis it is poissible to use ai(0) and 
a 2 (0) instead of a±(iu) and a 2 (iu), so (^) for r ^> A yields Casimir- Polder 
potential 



= _23 ai (0)a, W 
' r»A 4vr r 7 
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For A<jt. < r <C A (A ai . is atomic size) frequencies uj ~ 1/A give essential 
contribution to (B3). Then the last term in (|32|) gives main contribution and 



C/(r) 



do; ai{iu) a>2{ioj) . 



(34) 



Formulas (|33|) and (|34]) are non- dispersive and dispersive contributions to 
the energy of interaction of two neutral atoms. 



4 Energy of a dilute ball 

The change in the ground state energy of a dilute dielectric ball is generally 
equal 



dE = lv 



dE 

dV+ — 

e=const uE 



de , (35) 

V=const 



where V is the volume of the ball. Usually only the first term had been taken 
into account in the calculations of Casimir force and, as a consequence, the 
number of particles inside the ball wasn't conserved and the repulsive force re- 
sulted in most calculations. However, from physical considerations all atoms 
constituting the ball interact at large distances by attractive Casimir-Polder 
potential, so non- dispersive component of the force should be attractive. 

Conservation of the number of atoms inside the ball is equivalent to the 
conditions 

(e - l)V = const, de = - - ~^ — (36) 



We start from the first term contribution to fl35|) which we denote dE±. 
Figure 2 clarifies designations. When particles are inserted into the (infinitely 
small) shaded area near point B, they interact via Casimir-Polder (retarded) 
potential with particles inside a ball of radius a separated by distances r > A 
from the point B, i.e. they interact via Casimir-Polder potential with parti- 
cles inside the shaded region of the ball (regions J + II). Particles separated 
from the point B by distances X at . < r interatomic <C r < A interact with parti- 
cles near point B via van der Waals (non-retarded) potential. Obviously the 
situation is simplified for clarity here. Generally one has to take the concrete 
model of dispersion for dielectric permittivity, integrate interaction of par- 
ticles near point B via potential (|3l|) with particles separated by distances 
greater than interatomic distances from point B and calculate the energy 
change according to formula (|35|). As we are not intending to study specific 
model of dispersion here but rather want to get a suitable approximation for 
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Figure 2: contribution to dE\ 



non—disp. 



low-frequency contribution to the energy, we approximate area of integra- 
tion by its division on the retarded region (r > A) and non-retarded region 
{? 'inter ■atomic < t < A). According to discussion in the previous section, for 
distances r > A the frequency dispersion can be neglected, so this region 
in our approximation gives non-dispersive contribution to the energy , the 
region r interatomic < r < A gives main dispersive contribution to the energy. 

Area of integration for non-dispersive case is shaded on fig. 2, so in our 
approximation the non- dispersive contribution to dE\ can be written as 



dEi non-disp. = 47ra da J d r 

i+ll 



-23 a(r B , lo = 0) p(r B ) a(r, uj = 0) p(r) 



47r|rg — r 

\2 



7 



, o , (e(u = - 1 2 f -23 

= 4?ra 2 da v v - ' \ d 3 r — -, 37 

(47r) 2 J 47r|r s - r| 7 

/+// 

where p(r^) and p(r) are densities of particles in a unit volume. The integral 
over the region I can be simply obtained from ([21]) and arguments of the 
previous section. We just substitute a=^a — A, r=^a in formula (|2T| ) 
according to fig. 2 and obtain: 

-23 23 (q-A) 3 (5a 2 + (a-A) 2 ) 

4vr|r B -r| 7 15 a(2a - A) 4 A 4 K } 

i 



9 



Figure 3: contributions to dE 2 



non—disp. 



The integral (|37|) over the region II can be calculated analytically using (p0|): 

2-Kp 2 sin 9 



d 3 r = dp d6 — 

\ r B — r| y y (a 2 +p 2 — 2a»cosy)2 

II a-A 6»„ 



7 



2vrp V 1 1 tt(5A 4 + 192a 4 - 96a 3 A) 1 vr(5a - 4A) 
P ~ J P p6~ 480 c^A 4 30 a(2a - A) 4 

a— A A 

(39) 

From (|3"7|), ( |3BD and ( ]3"9"D after integration over a (keeping e = const) we 
obtain expression for Ei non _ disp : 

, ^o/ 23 7 23 5 23 1 \ 
E lnon . dlsp . = (e - 1) (-^^ + + 1536-J + const (40) 

Here V and S are volume and surface area of the ball respectively. Expression 
looks similar to formula ( fI5| ) obtained by G.Barton |TB[ in the framework 



of quantum mechanical perturbation theory. There is still a little difference 
between these formulas. At first sight it seems that in our approach it is 
impossible to find the term proportional to 1/ A and this term has no influence 
on the force. However, this is not so. The term proportional to 1/A is 
determined uniquely in our approach when taking into account the second 
term in (|35D and equation (|36|) . 

Now we turn to the second term in Q3"5D and calculate its non-dispersive 
contribution to the energy change denoted by dE2non-dis P .- Figure 3 illus- 
trates different areas of integration over variable r and analogous to fig. 2 in 
a sence - it is necessary to integrate Casimir-Polder potential between (in- 
finitesimally small area around) point G\ (C 2 , C 3 ) and shaded area of the 
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ball separated from point C\ (C2, C3) by distances r > A and then integrate 
over points C±, C2, C3 inside the ball. This term can be written as follows: 



dE. 



2 non—disp. 



' " ' ' ] f f dV(v c )dV(r) - , 23 - , (41) 



v=const An An J J An\rc — r\ 7 

\rn—r\>\ 



where and r are integrated over all points inside the ball of radius a 
satisfying the condition \rc — r| > A. 

To facilitate understanding of calculations in what follows we present 
major steps explicitly Let us denote 

g(r, p, p) = Anr 2 ■ 2np 2 ^— (42) 

rpp b 

and write dEmon-disp. (after comparison with (|39|) and fig. 3 everything 
should be clear): 

23e-lfe 
An An An 



x 



A A+r r+p A a r+p 



j dr J dp J dp + J dr J dp J dp^j g(r,p,p) see fig. 3.1 + 

A-r A r+A p-r 

a— A a p+r a— A r+A r+p 

+ I dr dp dp + dr dp dp + 



A r+A p—r 
a— A r— A r+p 



+ dr dp dp\g(r,p,p) see fig. 3. 2 + 



A r—p 

a a r+p a r— A r+p 



+ (^J dr j dp j dp + J dr J dp J dp^j g(r,p, p) see fig. 3. 3 



a— A r— A A a— A r—p 



23 V 23 S 23 1 23 1 \ 
= 26£{£ - 1 K-m^ + WAn^-WAnJ + \^n-a)- (43) 

The expression for non- dispersive force on a unit surface can be obtained 
using (|35|) and (j36|): 

1 (i-E 1 ^ 1 dE\ non —(iisp. dE^non—disp. _ 

tnon-disp. = -^2~-[^ = 7Z = 

23 / 4 A 1A 3 7 A 4 \ /JJX 
+ (44) 



1287T 2 A 4 V 3 a 2 a 3 48 a 4 
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This force is attractive for A < a , which reflects attractive nature of Casimir- 
Polder potential. 

To calculate energy, one has to take into account equation ( PS] ) and inte- 
grate force. As a result, 

^ 2 f 23 V _23_^_ _ 23 1 \ 

bnon-dis P . - (e I) { 128 7r 2 A 4 + 3 847r 2 A 3 3 84 vrA + 1536 vr J ' { } 

The term proportional to 1/A is determined uniquely now. Moreover, com- 
parison of ( |40D and fl45|) shows internal self-consistence of our approach. 



5 Conclusions 

The formulas for Casimir-Polder energy of interaction between a polariz- 
able particle and a dilute dielectric ball (|2T|),(|2"2"|) and non- dispersive contri- 
bution to Casimir surface force on a dilute dielectric ball (|4]) are derived 
analytically. Our approach directly shows correspondence of this force and 
Casimir-Polder potential. In section 3 this topic has been investigated from 
general viewpoint for arbitrary homogeneous dielectric media. The proof 
of equivalence in summation via Casimir-Polder - van der Waals type po- 
tentials and Casimir energy in the lowest (e — l) 2 order is the subject of 
section 3. This equivalence made it possible to calculate non-dispersive con- 
tribution to Casimir energy of a dilute dielectric ball by direct summation 
of Casimir-Polder potential. When the number of atoms constituting dilute 
ball is assumed constant during the ball expansion , the Casimir surface force 
is attractive and includes volume and surface terms in the order (e — l) 2 . 



6 Appendix 

For the benefit of the reader it is helpful to give here an overview of different 
approaches which have been used to derive the finite when A — > term 
Ef = 23/(15367ra) in the energy expression ( flop . 

First attempts to calculate Casimir energy of a dielectric ball gave a lot 
of different answers, as it is described in iTBl. The mathematical reason for 



some of these differences was found in the work |14[ , where by use of Debye 
expansion for Bessel functions and ^-function the correct limits on Ef were 



established. Later in the article the term Ef was calculated numerically 



with high accuracy, which made it possible to establish equivalence of Casimir 
effect and retarded van der Waals (or Casimir-Polder) energy for the case of 
non-dispersive dielectric ball. 
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Mutual retarded van der Waals energy for molecules inside the compact 



sphere (ball) was calculated by K.Milton and Y.Ng JT6| and its finite part 
after regularization in terms of gamma functions was first given by Ef in the 



work 16 



In the article [15| Casimir energy of a dilute dielectric ball was studied 
via the formula which has the following form for nondispersive case: 



E, 



c 



— 1 r f°° d 

Ly^(2Z + l) / dxx--F l 

57ra Jo dx 



(46) 



where 



F,(x) = — 



d 

A\dx 



(em, 



2 12 II 
X \S — SS 



12 II 

e — ee 



(47) 



and Riccati-Bessel functions are assumed to depend on argument x. 

The formula (f4*6|) can be studied using formula ( |T8| ) as well. We put 
k = l,p = 1 — m, where we are interested in the limit m — > 0. For Eq we 
find 



E, 



c 



lim 



(e 



If / 23 
V 1536 



(4f 



Recently the formula fl46|) has been studied in the mode summation 
method |T^], there divergent terms have been analysed and the finite re- 
sult Ef has also been derived by making use of addition theorem for Bessel 
functions. 

Another approach based on quantum mechanical perturbation theory was 



suggested by G.Barton in the work [[L3[]. The Casimir energy was obtained 
there in the form: 



E 



3 V 



+ 



2tt 2 A 4 
+ (^-l) 2 



V 



+ 



7 S 



1 1 



23 1 



128vr 2 A 4 360tt 3 A 3 20tt 2 A 15367ra' 



(49) 



Here V is the volume and S the surface area, 1/A is an exponential cutoff on 
wavenumbers. In this approach the contact terms haven't been subtracted, 
this is why the term proportional to e — 1 is present in (W). We have already 
discussed this formula in section 4, its structure is similar to fl45|) . The cutoff 
independent term is essentially the same in all approaches. 

The approach based on quantum statistical mechanics was developed in 
JIB]. This work led to similar results. 



13 



The theory of QED in a dielectric background was studied in || by path 
integral and (^-function methods, (^-function can be written as in ||: 

C ( s ) = ^r(21 + l)/ duu- 2s — ln(AA). (50) 

For a dilute ball it is possible to expand the logarithm in powers of e — 1, 
here we consider only the order (e — l) 2 of this expansion: 



C(s)= {£ / )2sin7rS f>/ + l) I" duu-^^ua), (51) 

where is defined in (f4"T|). Detailed derivation of the finite part of ( - 

function deserves a single article , here for completeness we only present the 
result for its finite part Q : 

C/M = -(«-!)■ ^Hlq" 2'-' ~ 3S + 4 >'j<- 23 + 2) , (52) 

with no poles in this expression for s < 1. From this expression we find 

cjtm = ^L^L, (53) 

; 2 1536 ira ( J 

we only note that even ^-function method meets severe problems with ana- 
lytic continuation and regularization in the case of a dilute dielectric ball. 

Most divergences and problems met so far in different approaches are 
related to short distance behaviour. In this article it is shown explicitly 
that such problems arise when parameter A — > 0. In dielectric media in the 
real world A ^ 0, this is why the main problem which still remains in other 
approaches is not to remove regularization, but to identify parameters of 
regularization in final expressions with physical quantities. 
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